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Abstract. 

We study the one-dimensional two-state totalistic probabilistic cellular automata 
(TPCA) having an absorbing state with long-range interactions, which can be 
considered as a natural extension of the Domany-Kinzel model. We establish the 
conditions for existence of a site-bond representation and self-dual property. Moreover 
we present an expression of a set-to-set connectedness between two sets, a matrix 
expression for a condition of the self-duality, and a convergence theorem for the TPCA. 
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1. Introduction 

Probabilistic versions of cellular automata can be considered as discrete-time Markov 
processes with parallel updating, which are useful in a large number of scientific areas 
[l]. As a special case, the Domany-Kinzel (DK) model introduced and studied by [2j is 
one of the simplest examples of the probabilistic cellular automata. The model is defined 
on a lattice, and the two different states for a site can be said as empty or occupied. The 
state of a given site depends on the number of the occupied states in its two nearest 
neighbors at the previous time step, where this property can be expressed with the 
totalistic rule |U EH E]. In the DK model, the transition probability from an empty 
neighborhood to an occupied state is zero, so the empty configuration is an absorbing 
state. The totalistic probabilistic cellular automaton (TPCA) is a natural extension of 
the DK model, where we extend the number of neighbors to a finite integer N(> 2). 
Remark that N = 2 case becomes the DK model. Concerning the DK model, various 
results are known (see our recent papers [El El IS] and references therein). As for N = 3 
case with two absorbing states, a rich phase diagram is reporeted by using mean-field 
approximations and numerical simulations (see Refs. [0100], for examples). However 
there are few rigorous results on the TPCA for a general N > 3. In this situation, the 
purpose of this paper is to give some rigorous results for a two-state iV-neighbor TPCA 
with an absorbing state. More precisely, we reveal a relation between the site-bond 
representation and the self-duality for the TPCA. Furthermore we present an expression 
of a set-to-set connectedness between two sets, a matrix expression for a condition of the 
self-duality, and a convergence theorem for the TPCA. It is known that the self-duality is 
a very useful technique in the study of stochastic interacting models. Because problems 
in uncountable state space (typically configurations of {0, 1} Z ) can be reformulated as 
problems in countable state space (typically finite subsets of Z). For some applications 
of the self-duality on discrete-time case (e.g., oriented bond percolation model) and 
continuous-time case (e.g., contact process), see Refs. [TTJ [TJJ EH Ej . So our results 
here might be useful in obtaining rigorous results on the phase diagram and critical 
properties of such models. 

We have organized the present paper as follows: Section 2 is devoted to the 
definition of our TPCA. In Section 3, we explain a site-bond representation and a 
self-duality for the TPCA. Moreover an expression of a set-to-set connectedness for the 
TPCA is given. Section 4 treats a matrix expression for the self-duality of the model. 
Section 5 contains a convergence theorem. Finally, concluding remarks and discussions 
are presented in Section 6. 

2. Definition of the TPCA 

First we give the definition of the DK model (A = 2 case). Let £^ C Z be the state of 
the process with parameters (^1,^2) £ [0, l] 2 at time n starting from a set A C 2Z. Its 
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P{x G = - 1) + + 1)), 

and given the events {x G Cn+i) are independent, where P(B\C) is the conditional 
probability that B occurs given that C occurs, and 

/(0) = 0, /(l)=j)i, /(2)=p 2 . 

If we write £(x,n) = 1 for x E and £(x, n) = otherwise, each realization of the 
process is identified with a configuration £ G {0, l} 5 with S = {s = (x, n) G Z x Z + : 
x + n = even }, where Z + = {0, 1,2, . . .}. The model with (pi,p 2 ) = (1,0) becomes 
Wolfram's rule 90. For more detailed information, see Section 5 in [T2*j . 

From now on, we introduce a long-range TPCA. In order to clarify the definition, 
we consider two cases of N = even and N = odd respectively, where N is the number 
of the neighborhood. 

(i) iV = 2L (L = 1, 2, . . .) case. The space of sites is denoted by So = {s = (x, n) G 
ZxZ + :x + n = even } and the space of bonds is 

B = {((x, n + l),(x-2L + l,n)),..., ((x, n + 1), (x - 1, n)), {{x, n + 1), 
(x + 1, n)), . . . , {{x, n + 1), (a; + 2L - 1, n)) : (a;, n + 1) G 5 }. 

For any initial set A C 2Z, the process £^ satisfies 

L 

k=-(L-l) 

and given the events {x G £^ +1 } are independent, where 

f(0)=p = 0J(l)= Pl ,...,f(2L)=p 2L , 

with pi,P2, ■ ■ ■ ,P2L G [0, 1]. This process is considered on the space So- The case L — 1 
is equivalent to the DK model. 

(ii) N — 2L + 1 (L = 1, 2, . . .) case. The space of sites is denoted by 

= i s = ( x i n ) G Z x Z + }, 
and the space of bonds is 

B l = {((x, n + 1), (x - L, n)), . . . , ((x, n + 1), (x - 1, n)), ((x, n + 1), 

(x, n)), ((x, n + 1), (x + 1, n)), . . . , ((x, n + 1), (x + L, n)) : (x, n + 1) G Si}. 

For any initial set A C Z, the process ^ satisfies 

L 

k=-L 

and given the events {x G £„ + i} are independent, where 

/(0) = = 0, /(l) = pi, . . . , f(2L + 1) = p 2L+1 , 
with Px, P 2, ■ ■ ■ ,P2L+i G [0, 1]. This process is considered on the space Si. 
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If p = < pi < p 2 < . . . < Pjv-i < Pn, then the A-neighbor TPCA is said to be 
attractive. If not, it is said to be non- attractive. The attractiveness means that having 
more particles at one time implies there will be more particles at the next time(s). 

In general, few rigorous results on the non-attractive model are known compared 
with the attractive model. Much more informations on the results for the non-attractive 
N = 2 case are shown for instance in [7 r So studying a general model including the 
non-attractive case, such as is the one in this paper, is very important. 

3. Site-bond representation and self-duality 

From now on we consider only (pi,p2) G D* case, where the subset of parameter space 
.D* is defined in the following: 

D* = {(pi,p 2 ) : < p x < 1, < p 2 < 1 andp 2 < 2pi}. 

A reason for introducing the set is that for any point outside of D*, it will be shown 
that there exists no site-bond representation for the TPCA (see Proposition 1 below). 
Suppose that pi,p2 G D* are given. We put a = p\j{2p\ — p 2 ), (3 — 2 — p 2 /pi. Let 
(S,B) = (Sq,B ) or (Si,Bi) and R denote the set of real numbers. If there exist 
a, (3 G R such that p n G [0, 1] is given by 

p n = a[l-(l-pn (1) 
for n = 3, 4, . . . , N, then we call "an iV-neighbor TPCA with {p n : < n < N} has a 
site-bond representation with a and In fact, when a, (3 G (0, 1], the TPCA can be 
considered as an oriented mixed site-bond percolation with a long range interaction in 
the following way. On the space (S,B) we define 

X(S) = {0,1} S , X(B) = {0,1} B , X = X(S)xX(B). 

For given ( = ((1,(2) G X, we say that s = (y,n + k) E S can be reached from 
s' = (x, n) E S and write s' — > s, if there exists a sequence s , s 1; s 2) ■ ■ ■ ; s k of members 
of S 1 such that s' = so, s = s^, and 

Ci(sj) = l,z = 0, l,...,fc; ( 2 ((s h s i+1 )) = l,i = 0, 1. 

we also say that G C S can be reached from G' C S and write G' — ► G, if there exist 
seG and s' G G' such that s' — > s. 

We introduce the signed measure m on X defined by 

m(A) = a kl {l - a) n f3 k2 (l - f3) n , 

for any cylinder set 

A = {(Ci, C2) G X: Ci(si) = l,i = l,2,...,h; Ci(sJ) = 0, 2 = 1, 2, . . . , Ji; 
C 2 (6,) = l,z = l,2,...,fc 2 ; Ca(6J) = 0,z = l,2,...,j 2 }, 

where s 1; s 2 , . . . , s^, s 2 , . . . , s'j are distinct elements of S, and 61, 6 2r . . , b k ,b[, b' 2 ,. . . , b j 
are distinct elements of B and a = p\j (2pi — p 2 ), [3 = 2— Pijp\- 
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We should remark the next three facts. If p 2 < 2p\ and p 2 > 2p\ — p\, then a > 1 
and (3 G (0, 1). If p 2 < 2p^—p\ andp 2 > pi, then a,j3 G (0, 1]. Moreover, if p 2 < p 1 , then 
a G (0, 1) and (3 G (1,2]. From the above observation, we see that the measure is not 
probability measure in the first and third cases, since a > 1 and 1 — (3 < respectively. 
See Refs. |S] and [Zj for more details on the site-bond representation. 

It is noted that if an iV-neighbor TPCA has a site-bond representation with a and 
(3, then a = p\J {2pi — p 2 ), [3 = 2 — p 2 /p\, since Eq. ((H) gives 

P! = aP, p 2 = a[l-(l-f3) 2 }. 

As we will see in the proof of Proposition 1, the following result can be obtained: an 
iV-neighbor TPCA having a site-bond representation with a and (3 is attractive if and 
only if pi < p 2 . In other words, an iV-neighbor TPCA having a site-bond representation 
with a and (3 is non- attractive if and only if pi > p 2 . Then we have 

Proposition 1 We consider iV-neighbor TPCA with {p n : < n < N}. Assume 
{P11P2) G .D*. Let a = pj/(2p 1 — p 2 ) and (3 = (2pi —p 2 )/p\. (i) pi > p 2 (non-attractive 
case) andp n = a[l-(l-(3) n ] (0 <n<N), thenp n G [0,1] (0 <n<N). (ii) p x < p 2 (at- 
tractive case), p n = a[l-(l-/3) n ](0 < n < N), andp N < 1, thenp n G [0, 1](0 <n<N). 

Proof. First we consider a relation among p = 0,p±, . . . ,Pn- By using p n = 
a[l - (1 - (3) n ) for n = 0, 1, . . . , JV, we see 

p n -p n ^ = a(3{l-(3) n -\ (2) 
Pn - Pn _ 2 = a{l - f3) n ~ 2 {l - (1 - (3) 2 }. (3) 

Recall that 

1- (3 = (p 2 -pi)/pi, a(3 = p 1 . (4) 
Now we consider > p 2 (non-attractive) case. In this case, Eq. (0J) gives 

-1<1-P<0, 0<ap<l, 0<ot<l. (5) 
From Eqs. (J2J), ©, and ©, we see 

p = < p 2 < PA < P6 < ■ ■ ■ < P2n < P2n-1 < ■ ■ ■ < Ph < P3 < Pi < 1 • 

Thus we obtain p n G [0, 1] for any n = 0, 1, . . . , N. Next we consider p\ < p 2 (attractive) 
case. It is easily checked that p n < p n +\ for any n = 0,1, . . . , N — lby using Eq. (J2J) 
and < 1 — (3 < 1. So we have 

p = < p 1 < p 2 < p 3 < . . . < p N -i < p N . 

Therefore a necessary and sufficient condition for the site-bond representation in this 
case is p^ = a[l — (1 — f3) N ] < 1. The proof of Proposition 1 is complete. 

Remark that p n = a[l — (1 — f3) n ] — > a as n — > 00, since |1 — /3\ < 1. Noting that a < 1 
if and only if p 2 < 2p 1 — pi, we conclude that a sufficient (and not necessary) condition 
for a site-bond representation of an iV-neighbor TPCA is p 2 < 2p\ — p\. Furthermore 
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Corollary 2 If an iV-neighbor TPCA with {p n : < n < N} and iV > 3 has the 
site-bond representation, then p n can be expanded with pi and p 2 as follows: 

(P2-P1T' 1 (P2-Pl) n ~ 2 . , (P2~Pl) 2 , 
Pn = ;j=2 + ^3 + • • • + + P2, 

Pi Pi Pi 

for any n = 3,4, . . . , N. 

Next we define a set-to-set connectedness for the TPCA from a set A to a set P by 
a(A,B)= lim P(£fnP^0), 

+00 

if the right-hand side exists. As in a similar argument of N = 2 case (the DK model) 
jS], we can easily extend the case to a general N- neighbor TPCA: 

Proposition 3 Let (pi,P2) £ P*- We assume that an iV-neighbor TPCA with 
{Pn '■ < n < N} has a site-bond representation with a and /3, where a = p\j (2pi — p 2 ) 
and /3 = (2p\ — p 2 )/pi. Then for any A with \A\ < 00 , we have 

a(2Z,A) = a\ D \l-a)\ A \ D \a{D,2Z), 
when N = 2L, and 

when AT = 2L + 1. 

From now on, we consider a self-duality of the TPCA. An iV-neighbor TPCA £ n is 
said to be self-dual with a self-duality parameter x if 



£ ( x \tfnB\\ E f x \(»nA\\ ( n = ,l, 



holds for any A, P C Z with |A| < 00 or |P| < 00. The above equation is called a 
self-duality equation. Then we have 

Theorem 4 Let (pi,p 2 ) G P* with a[l - (1 - /3) N ] < 1. An A^-neighbor TPCA 
with {p n : < n < A^} is self-dual with a self-duality parameter (a — l)/a if and only 
if this model has a site-bond representation with a and (3, where a = p\j{2p\ — p 2 ) and 
(3 = (2pi-p 2 )/px. 

Proof. Theorem 1 in implies that an A^-neighbor TPCA with transition prob- 
abilities {p n : < n < N} is self-dual with a self-duality parameter x is equivalent 
to 

(piX + q { y = {pjX + qjY, (6) 
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for < i, j < N, where = 1 — pj. It can be confirmed that x = (a — l)/a with 
a — Pi/ (2pi — P2) satisfies Eq. (0). Moreover we see that Eq. (0 if and only if 

(pix + qi) n = p n x + q n , (7) 

for < n < N. In fact if we take i = 1 and j = n, Eq. (0 becomes Eq. (J7J). Conversely 
we have 

(pix + = (pix + gi) 4 -' = (pjX + qj) 1 . 
Therefore we put y — x — 1. Then Eq. (JJJ) becomes p n y + 1 = (pi?/ + l) n . So we see that 



r=l 



E 



r=l 



r-l 



p 2 - ggi \ 1 A AA / p 2 ~ 2gi x ' 
Pi J t^i W V Pi 



E 

r=0 



n> 1-^ 

77 V Pi 



= a[l - (1 + (pa - 2p 1 )/p 1 ) n ] 
= a[l-{l-0) n ]. 
So the proof of Theorem 4 is complete. 

From this theorem it can be seen that two conditions for the self-duality and for 
the site-bond representation are equivalent under an assumption that (^1,^2) £ -D* with 
p JV = «[l-(l-/3) Ar ] < 1. 

4. Matrix expression 



In this section we consider a criterion for the self-duality based on a matrix expression. 
Let 



and 



X(x) 



1 1 

1 x 



X N (x) = X(x 



,®iV 



where <S> indicates the tensor product. For example, in the case of N = 2, we have 



X 2 (x) = X(x)®X(x) 



1111 

1 X 1 X 
1 1 X X 

1 X X x 2 
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For fixed i,j G 1,2, .. . ,2^, the values of ii,i 2 , • • . , ijv and ji, J2, ■ ■ - iJn £ {0, 1} are 
defined by the binary expansion of « and j in the following: 

i — 1 = %\l + 222 + h«7V-i2 +«at2, 

J-l=Ji2 +J 2 2 H hjAr-i2 + j N 2 , 

and / = «i + ?2 + • • • + ijv- Moreover we introduce 2 W x 2 W matrix Pjv whose 
element is defined by 



Jh) (32) 
Pi Pi 

where jfp = pj, pf^ = 



. . . T) (ijv) 

Pi > 

1 — pi — qr. For N = 2 case, the above definition gives 



1 











Qi PiQi PiQi pi 
Qi P1Q1 P1Q1 Pi 

ql P2<?2 P2<?2 P\ 



since p — 0, q = 1. Then we have the next result: 



Theorem 5 A necessary and sufficient condition of the self-duality for an iV-neighbor 
TPCA with {p n : < n < N} with a self-duality parameter x is that x satisfies 

P N X N (x) = '(PnX^x)), (8) 

where t means the transposed operater. 

Proof. Let c^\x) denote (i,k) element of PnXn(x). First ii, i 2 , ■ ■ ■ , in and 
k u k 2 ,...,k N are defined by i - 1 = i^' 1 + i 2 2 N ~ 2 + ••• + i N ~i2 l + i N 2° and 

k-1 = k 1 2 N - 1 + k 2 2 N ~ 2 -\ hfcjv_i2 1 + A;j V 2 , respectively. We put I = h + i 2 -\ Vi N 

and K = ki + k 2 + ■ ■ ■ + /cat. Moreover {k\, k 2 , . . . , k^} is divided into two sets such as 
{k mi ,k m2 , . . . k mK :k u = l for any u E {mi,m 2 , . . . ,m K }} and {k h ,k h , . . . k lN _ K : k v = 
for any v £ l 2) . . . , In-k}}- Then we see that 



x) 



G{0,1} G{0,1} 



Pi 



• Pj x Pi 



x x 



Pi 



•Pj X 



= £ 

im 1 ,.-.,im if €{0,1} 

/ (1) , (0) \K 

= (Pix + qif. 
where x^ = x and X( ) = 1. In a similar way, we have 
c kP( x ) = (PrX + Qk) 1 - 

So we have 

P N X N (x) = '(PnX^x)), 



(jmj ) 



X 



(jm K ) 



X X f - 



Oh) 



■ X 



til 
■Pi 

k lN- 
til N 



•) 



■ X 



tim K ) 
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if and only if 

for 1 < i, k < 2 N . Forthermore it is shown that the last equation is equivalent to 

(jo/x + qi) K = (p K x + q K y, 

for 1 < J, K < N. Combaining the above result with Theorem 1 in [15J, we obtain the 
desired conclusion. 

For N = 2 case (the DK model), Theorem 5 gives 
(Pix + qi) 2 = p 2 x + q 2 . 
This result was shown in E3 CHI by using different methods. 

5. Convergence theorem 

In order to obtain a convergence theorem for the TPCA with {p n : < n < N} having a 
site-bond representation with a = p\j{2p\— p 2 ) and j3 = (2pi—p 2 )/pi, which corresponds 
to the result given by 0, a new process r\ n is introduced as follows. For simplicity, we 
first consider the DK model (N = 2 case) with (pi,p2) € -D*. Put p* = max{pi,p 2 }. A 
new process defined below is called p^-DK-dual. We can see the thinning-relationship 
by coupling the DK model and p*-DK-dual. We split both the DK model and the p*- 
DK-dual into two phases, and we will allow the first phase to occur at times n + (1/2) 
for n G Z + . 

(i) Let n be the distribution of the p*-DK-dual at time 0. 

(ii) At time n — 1/2, it undergoes a p*-thinning. In general, for p e [0, 1], the p-thinning 
of a set A C Z is the random subset of A obtained by independently removing each 
element of A with probability 1 — p. 

(iii) Start the DK model at time n = with the same configuration as the p*-DK-model 
(which is defined by /(0) = 0, /(l) = p\/p*, f{2) = p 2 /p*) at time n = 1/2. 

(iv) Couple the processes together until time no — (1/2) for the DK model, n for the 
p^-DK-model. This can be done because the transitions for the DK model are the 
same as those for the p*-DK- model lagged by time unit 1/2. 

(v) Perform a p*-thinning for the DK model at time no- 

The distribution of the DK model started and ended as a p^-thinning of the p^-DK- 
dual. As in the DK model, we can define a new process r] n as p^-TPCA-dual for the 
A^-neighbor TPCA, where p* = m&x{pi,p 2 , . . . ,Pn}- Recall that the following duality 
holds for the TPCA £j™ and the p*-TPCA-dual rj^ starting from A, (see Theorem 2 (2) 
in H5|). 
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Theorem 6 (|15j) We assume an iV-neighbor TPCA with {p n : < n < N} hav- 
ing a site-bond representation with a = p\/(2p 1 — p 2 ) and (3 = {2pi — p 2 )/pi- Suppose 
{PiiP2) £ D*. For any A, B with \A\ < 00 or \B\ < 00, we have 



for any n > 0, if x = 1 — (2pi — P2)p*lp\- 

Note that if pi > p 2 (non-attractive case for the TPCA with a site-bond represen- 
tation), then p* = p\. By using the same argument for the DK model in Section 7 of 
jH], we can get the following convergence theorem for the TPCA: 

Theorem 7 Assume p\ > p 2 with (^1,^2) £ D*. If the initial measure v of the TPCA 
is a.s. (almost sure) infinite, then we have 



as n — >• 00, where the limit measure is uniquely determined by E((— (p\— p 2 ) / Pi)^ 71 ^) = 
P(\v^c\ — 0) for any A with \A\ < 00. 

6. Conclusions and discussions 

In this work we have presented rigorous results on the site-bond representation, the 
set-to-set connectedness, the self-duality, the matrix expression, and the convergence 
theorem for an iV-neighbor TPCA with {p n : < n < N}, where po = and N > 2. 
An interesting feature of our model is that the dominant parameters on some properties 
are only pi and p 2 among {p n : < n < N}, see Proposition 1 and Corollary 2, for 
examples. 

Arrowsmith and Essam J7| gave an expansion formula for a point-to-point 
connectedness for the oriented mixed site-bond percolation, in which each term is 
characterized by a graph. Konno and Katori jfj] extended this formula to iV = 2 case 
(the DK model). So it is shown that site-bond representation, self-duality, and the above 
graphical expansion formula hold in the DK model. Thus one of the interesting future 
problems is to extend the relation to a general iV-neighbor TPCA considered here. 

Finally we mention a relation between our discrete-time model and a continous- 
time one corresponding to it. In the continuous-time case, an infinitesimal generator for 
the model corresponds to Pn in Theorem 5. In fact, Eq. (JHJ) in Theorem 5 corresponds 
to Eq. (17) in JH] for N = 2 case (an extension of their result to an iV-neighbor case 
can be easily obtained). Recent works on the duality for continuous-time interacting 
particle systems are presented in (THl E01 EH E2J • 
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